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The downside risk of a portfolio of (equity) assets is generally substantially 
higher than the downside risk of its components. In particular in times of crises 
when assets tend to have high correlation, the understanding of this difference 
can be crucial in managing systemic risk of a portfolio. In this paper we gen- 
eralize Merton's option formula in the presence jumps to the multi-asset case. 
It is shown how common jumps across assets provide an intuitive and powerful 
tool to describe systemic risk that is consistent with data. The methodology 
provides a new way to mark and risk-manage systemic risk of portfolios in a 
systematic way. 



Introduction 



It has been argued that one of the factors that triggered the downfall of Long 
Term Capital Management (LTCM) was its failure to incorporate "fat tails" 
of asset price distributions properly into investment decisions as well as risk 
management [1]. Today financial institutions systematically deduce fat tails 
from the option markets and incorporate this information consistently into the 
pricing as well as the risk management framework. Despite this progress it is 
interesting to note that little or even no efforts have been made in the marking 
and risk-managing of "fat tails" when individual assets are combined into a large 
portfolio. However a large part of a portfolio's fat tail is not the result of the 
fat tails of the individual components. It is this difference that can be associ- 
ated with systemic risk of a portfolio that will generally dominate investment 
decisions in times of crises. 

As pointed out [2] only roughly 50% of the downside risk (fat tail) of a portfolio 
of stocks such as the DAX, is due to the fat tails of its members. In a sense by 
ignoring the additional 50% of downside, financial markets have not yet fully 
learned the lessons from the LTCM debacle. In fact the current financial crisis is 
haunting us with the consequences of underestimating systemic risk in portfolio 
management. 

In this paper we argue that financial institutions should extend their market 
data to account for systemic risk not only in the area of credit but across all 
assets classes whenever assets are combined into a portfolio. The proposed "ex- 
tension of market data" would be similar, in some respect, to the generalization 
of at-the-money volatilities to all strikes. 

The first objective of this paper is to show that the steepness of the DAX skew 
compared to its components can be linked to the systemic risk of the "DAX 
portfolio". This is important because it allows us to infer the systemic risk 
from the market consistently. The second objective of this paper is to explicitly 
construct a parameterization that describes the systemic risk of large equity 
portfolios. The latter can serve as a candidate for the extension of the market 
data mentioned above. 

Recently, efforts have been made to extend Dupire's local volatility framework 
to incorporate state-dependent correlations into the dynamics of a portfolio [2, 
3]. While this leads to a numerically efficient generalization of Dupire's model, 
there are two main drawbacks. 

First, the model requires an a priori knowledge of the basket skew as an in- 
put. This is however not known in many cases. The method presented in this 
paper provides a natural way to interpolate skews to sub-baskets as well as 
cross-index baskets. 
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Second, while the extension of Dupire's model to local correlation consistently 
describes the correlation skew, it does not attempt to "explain" its causes. Con- 
sequently, it cannot be viewed as a straightforward methodology to monitor 
systemic risk. The approach presented in this paper defines a copula that ex- 
plicitly links the correlation skew to systemic risk and hence could be a can- 
didate to improve the (tail) risk-management of large portfolios. Work in this 
direction has been pursued by others [4], [8]. The paper is organized as follows: 
The first section reviews Merton's option formula that describes options in the 
presence of jumps. We then generalize this formula to the multi-asset case and 
show qualitatively how correlation skew can be induced in this framework. The 
next chapter generalizes the single-asset Merton formula in order to properly 
reproduce market prices of options at all strike prices. Finally we put both gen- 
eralizations together and define the "Merton-Copula basket" that enables one 
to define systemic portfolio risk in a proper way. 



A Merton basket option formula 



In 1976, Merton generalized the Black-Scholes formula to situations where jumps 
are present [5]. In this case the dynamics of the stock price S t is given by 



where r t ,qt,\dn\,a,w t denote the short rate, dividend yield, jump intensity, 
jump measure, diffusive vol and a Brownian motion respectively. The jump size 
Y — 1 has a mean k and a lognormal volatility +S, i.e. 



where z e TV (0,1). Equation 1 states that the dynamics of St is log-normal 
with occasional jumps occurring at rate A and with size Y — 1. 

As the jumps are assumed proportional (rather than additive shifts) as well 
as fully uncorrelated to the diffusive part of the dynamics, the Merton price 
of a call option Call(T,K) with maturity T and strike K is simply given by an 
infinite sum of Black Scholes prices, i.e. 



St 



= {r t -q t - \E X [Y - 1]) dt+(Y- 1) dn x + adwt 



(1) 




Call (T, K) = J2 



n=0 
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where 



BS(F,T,K,a)=Df(FN(d 1 )-KN(d 2 )) 



di/2 = 
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denotes the standard Black-Scholes option formula. 



It is very easy to interpret Eq.2 in intuitive terms: e ^ / is the proba- 
bility of n jumps happening until maturity. In this case the spot has suffered an 

extra depreciation of ^1 + kj on average due to the jumps which is reflected 
in the adjustment to the forward price in the Black-Scholes formula. 

The Merton model describes a simple way to incorporate volatility skew into the 
dynamics without loss of analytic tractability. However the model is too simple 
to give an exact fit to the observed volatility skew in practice. In Fig.l we plot 
the implied volatility skew for different levels of the diffusive volatility a. As a 
increases the jump generates less skew which is consistent with intuition. 
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Fig.l: The volatility skew for different levels of diffusive vol. The effect of 
jumps on the skew becomes less and less pronounced for high values of the dif- 
fusive vol. 

In the Appendix 1 we show that a way to retain a volatility skew for large 
a stocks is to scale the jump size with the diffusive vol itself, e.g. 




(4) 



i7oj k denotes a "vol-scale" and "vol-scale elasticity" respectively. 
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In order to generalize Merton's formula to basket options, Eq.l can be gen- 
eralized to N assets that evolve according to 



dS^_ 
5« 



= (r t - qi ] - XE X [y« - l] dn A ) + aWdwP + (V W - l) dn Xl 



dwf\dw[ ]) ) = p% ffusive dt i = l,...,N (5) 

Note that, despite the multi-assets character of Eq.5 only a single Poisson pro- 
cess is "employed" just like in the case of equation 1. As we show later, this fact 
will prove crucial in order to generate correlation skew as well systemic risk of 
a portfolio. 

Eq.5 states that the dynamics between jump events is governed by a multi- 
asset Black-Scholes dynamics with a diffusive correlation p^ff uslve 

In the following we define a basket 

n 

B t =J2 a i S t ] (6) 

i=l 

with fixed basket weights a%. If Q denotes the pricing measure, Df the discount 
factor and fcj the jump size of asset i, then a call option on a basket can be 
calculated according to 

{Basket) - Call (T, K) = Df E Q {{B t - (7) 
£ e- XT{ ^TMPricer (f\ (l + h)" e^ T , ...,F N (l + fc^)" e~^ T ,pf* s 



n=0 



Note that contingent to the number of jumps, the dynamics is log-normal and 
can hence well be approximated by a well known Black-Scholes three-moment 
basket pricer denoted by TMPricer (F\, . . . , Fjy 7 p)( see for example [6] ). 

Equation 7 can be viewed as the generalization of Merton's formula to options 
on a basket of assets. 

In Eq.5, 7 all assets "share" the same Poisson process defined by 

(A,ft,<f) • This 

triple defines a Poisson process with intensity A, a universal jump size k and 
jump volatility S. 

Individual jump-sizes ki are computed from the "universal jump-size" by means 
of 



In the following we present results of Eq.7 in combination with Eq.8 when ap- 
plied to the components of the DAX. Fig. 2 shows the Merton basket implied 
vols that are inferred from Eq.7 versus the vol of the DAX for ^A, k,S, Co 

(25%, -16%, 18%, 18%) 
Table 1 shows how parameters change with the time horizon under considera- 
tion. 
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Fig. 2: Reconstruction of DAX skew from its components: Merton basket 
implied vols versus market for (A, jfe, 5, a ) = (25%, -16%, 18%, 18%) 




Maturity A 


k Sk 


O"0 


1 Year 0.25 

2 Year 0.25 

3 Year 0.25 


-0.13 18% 
-0.137 15% 
-0.13 18% 


0.12 
0.12 
0.127 



Table 1 : Jump parameter dependence as a function of time that is required 
to reconcile the skew of the DAX from its components. The parameters appear 
high relative to historic values. However they, in parts, reflect risk-premia that 
are associated with the sale of deep out-of-the-money put-options. 

It is remarkable that the simple formulas of Eq.7, 8 leads to quite good agree- 
ment with the skew of the DAX. We remind the reader that typically only 50% 
percent of the DAX skew can be explained in terms of the skew of the underlying 
assets [2]. Hence the idea of "sharing" a Poisson process across different assets 
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seems an idea worth investigating as a means to create copulas that are consis- 
tent with data. The intuition behind this concept is easily explained: When a 
far out-of-the money put on a basket ends up in the money at time of maturity, 
it is likely that at least one jump has occurred during the life of the option. 
As jumps occur simultaneously across all assets in our model, they effectively 
increase the implied correlation for out-of-the-money puts giving rise to corre- 
lation skew. The situation is reversed for out-of-the money call options that are 
likely to pay out only when no jump event has occurred. The implied correla- 
tion is expected to be close to the diffusive correlation in this case. For 6 = 
the relationship between the diffusive correlation p^ff uslve an( j total expected 
correlation pij can be calculated explicitly with the result 



Note that in the limit of very frequent jumps where A — > oo the correlation goes 
to one confirming our intuition. Eq. 9 plays a useful role when calibrating the 



constant in order not to change the levels of the at-the-money basket vols for 
different choices of A during calibration. Higher levels for the intensity require 
lower diffusive correlations which can be deduced from Eq.9. 

Even though this explanation is highly plausible, it is not the only mechanism 
that yields tp correlation skew. For example, multi-asset stochastic volatility 
models that exhibit negative cross-correlations between spot returns of one asset 
with the vol of another asset generate cross- vanna sensitivity and hence correla- 
tions skew. It would be interesting to combine this effect with the one outlined 
above to a more complex model. The latter would decrease the intensities as well 
as the jumps-sizes of Table 1 that are required to achieve consistency between 
the index skew and the skew of its components. In many respects such a model 
could be considered the next step in the devlopment that is somewhat similar to 
progress in the modelling of FX derivatives where two different mechanism for 
the volatility skew- local vol and stochastic vol get combined to a more realistic 
description of this smarket. 

We summarize this chapter by stating that the ideas presented seem worth- 
while perusing. However the Merton models as well as the Merton basket for- 
mula need to be extended in order to make the calibration to the individual 
stock skews exact. 



dif fusive 
Pij 






needs to be kept 
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Single asset Merton extension 



The main drawback of Merton's model is its failure to calibrate option prices 
exactly to the market. An example for the mis-calibration is given in Fig. 3. 
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Fig.3:Example for mis-calibration of a single stock skew in the Merton model. 

In order to fix this problem we generalize the Merton Single asset formula of 
Eq.2 to 

Call (T, K) = 

£ e"- { -^f BS (f (l + Xk) H e--, T, K, f^f^j (10) 

The only difference is that the diffusive vol itself becomes a function of the 
strike price. The idea is that for a given set \X, k, S^j the calibration mismatch 
is picked up by a (K). 

The function a (K) can easily be determined by means of the fixed point algo- 
rithm [6]: The The algorithm starts with a constant zeroth guess a for the dif- 
fusive volsurface (Kj) = a for M different values of strike Kj, j = 1, . . . , M. 
Iteration i + 1 is obtained from iteration i according to cr^ +1 ) (Kj) — (Kj) + 
Act,- where A<jj is the vol-mismatch between the implied vol inferred from Eq.10 
using (Kj) and the market implied vol. For reasonable first guesses the 
algorithm typically converges very fast as the following diagrams show 
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Spotbased Strike Diffusive Vol Market Vol Mispricing 
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F&Ed Point algorithm 1 iteration results of diffusive vol of DBK&i 
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Fig. 4: First iteration results for the fixed point algorithm: 
Second iteration results for the fixed point algorithm: 



Spotbased Strike Diffusive Vol MarketVol Mispricing 

0.51 40.9% 52.6% 0.1% 

0.63 41.4% 48.4% 0.1% 

0.79 40.4% 44.2% 0.1% 

0.99 38.7% 40.6% 0.1% 

1.23 37.1% 38.4% 0.1% 

1.54 37.3% 38.2% 0.0% 
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Fixed Point algorithm 2 iteration results of diffusive vol of DBKQi 
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Fig. 5: Second iteration results for the fixed point algorithm: 

Equation 10 together with the fixed-point algorithm for the determination of the 
strike dependent diffusive skew defines the our "Merton single asset extension 
formula" for the pricing of options in the presents of jumps as well as "diffusive" 
contributions to the skew. 

It is important to realize that one can easily construct a stochastic process 
explicitly that yields Eq.10 in the presence of "diffusive skew" u(K) for exam- 
ple: 

- E i = (r t -q t - XE X [Y - 1]) dt + (Y — 1) dn x + a {X u t) dw t (11) 

where a (X t ,t) can be thought of a some kind of "local vol" with the difference 
that it is not a function of the spot St but rather the variable Xt = Sk where 
n is the number of jumps that have occurred up to time t. 
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The Merton Copula Basket 



The goal of this chapter is to incorporate the results of the previous chap- 
ter into a basket options formula along the lines of Eq.7. However this is 
slightly more difficult as the three-moment approximation that was applied in 
Eq.7 is not appropriate in this case. This is because contingent to a specified 
number of jumps the dynamics is not log-normal anymore but given by a lo- 
cal volatility of Eq.ll. In the following we outline the generalization of the 
three-moment pricer that was implicit in Eq.7 to "infinite" moments denoted by 
BasketCopulaPricer (Si,K, p, cn{K)) 

• Step 1: Dial a set of correlated normal deviates Wi i — 1, . . . , TV from a 
correlation matrix pij and convert them into uniform deviates according 
to Uj — N (vjj ) 

• Step 2: Invert the individual asset (option implied) distributions accord- 
ing to 

to obtain a set (ui) , (u 2 ) , ■ ■ ■ , (w/v)) 

• Step 3: Calculate payoff = (Jt a i S t ] ~ 

• Step 4: Go to step 1 Npath times, then average results and obtain 

Call (T, K) = BasketCopulaPricer (Si, K, p, <J l {K)) 

We are now in a position to define the Merton Copula basket price in the pres- 
ence of jumps as well as "local volatility" (see Eq.ll) 

Call(T,K) = J2 e ~ X1 BasketCopulaPricer (St (l + e^ T , K, pfj ffu * ve , a t (K)^ 

(12) 

where cri(K) is given by the solution of the fixed point algorithm. 

Fig. 6 shows the results in the case of A = recovering the regular Gaussian 
Copula result. As mentioned in the beginning of the paper that not even 50% 
percent of the skew of the basket is recovered in this case from the skew of the 
individual assets. 
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Fig. 6: The Merton Copula basket reduces to the standard Gaussian Copula 
in the case of A = 0. Less than 50% of the skew of the basket is recovered in 
this case. 

Fig. 7 shows the results after calibration. Note that for all choices of the calibra- 
tion set ^A, k, S, nj the approach is consistent with the skews of the individual 
assets. 
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Fig. 7: Calibration results of the DAX skew from its components for the set 
'A, k, 5, k) = (0.35, -0.16, 18%, 1, 17, 18%) . 
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Hence the set ^A, k, S, kJ describes "orthogonal new market-data" that quan- 
tifies the correlations skew. One expects a different set for each maturity. 



13 



Summary 



We presented a generalization of the Merton Jump Option formula to multi- 
assets as well as diffusive skews and found that common jumps are an appro- 
priate and intuitive way to define copulas that describe the systemic risk of a 
portfolio. Financial institutions should start treating systemic risk as "market 
data" across asset classes and start to mark, risk-manage and price these effects 
in a systematic way. For each maturity a tuple ^A, k, S, hij can serve as a candi- 
date to describe this systemic risk in an appropriate way. It is straightforward 
to generalize our framework to more than one Poisson jump in the case different 
sources of systemic risk are encountered. This makes the framework well suited 
to mark skews of sub-baskets as well as cross-index baskets. 



Appendix 1 



In this Appendix we motivate the scaling of the jumpsize by the volatility a 
according to Eq.4 for k = 1 in order to retain the volatility skew for large a. 

From Eq.2 the price of a Digital-Put option of strike K and maturity T 
Digital - Put(K,T) = E 13 (1 St <k) is given by 



OO 



Digital - Put (T, K) = ^ e~ AT ^— ^-7V(-d 2 ) (13) 
Hence the price change of the digital with a change in intensity A is given by 



dDigital- Put(T,K) ^ _ XT {\T) n Zn- F(] ' 



£e-^ M2 ,^=^ (14) 



n=0 v 

where ip(x) denotes the normal density.This is due to the fact that for small 
values for A the relationship between implied and diffusive vol can be approxi- 
mated by a 2 (implied) — a 2 (dif fusive)+Xk 2 and the last term in the numerator 
of Eq.3 is generally small for suffently out-of-the money options. Hence if one 
wants to retain the skew for large values of a, the factor in Eq.14 suggests to 
scale the jump-size according to fc— > ok. 
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